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■ We calculate the first correction beyond proximity force approximation 

, for a cylindrical graphene sheet in interaction with a flat graphene sheet or a 

Qh! dielectric half space. 

c : 

1 Introduction 

FT 

! The interaction of material bodies at distances in the micrometer and nanometer 

scale is of significant actual interest in view of applications in nanotechnology and for 
precision measurements of the Casimir force. The basic method for the calculation 
of the interaction forces is the Lifshitz formula pQ. It holds for planar or stratified 
material bodies having permittivity. Permeability as well as frequency dispersion 
and finite temperature may be included too. The generalization to non planar 
geometry is more complicated in case the variables in the underlying wave equation 
do not separate and for the interaction of two material bodies one is bound to the 
case of plane parallel geometry. However, most force measurements are done with 
a sphere or lens in front of a plane and there is a demand for generalizations to non 
planar geometry. 

The most frequently used method for non planar geometry is the proximity force 
approximation (PEA). Here the Casimir or van der Waals force known from planar 
geometry is taken at the local distance and then averaged over the surfaces. In this 
method the curvature of the surfaces in included in the sense of a first correction. 
The method holds obviously for small deviation from the plane parallel geometry. 
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Its precision cannot be estimated since it is impossible to calculate higher order 
corrections by this method. 

The method of PFA dates back to Derjaguin [2] and only recently it became 
possible to go beyond analytically. In [3] the first correction beyond PFA was cal- 
culated for a cylinder in front of a plane with conductor boundary conditions. For 
the part containing Dirichlet boundary conditions this was confirmed numerically 
by the world line method [I] . The opposite case of large and medium separations is 
easier and more detailed results are available, see [5] and [6]. 

Another line of generalizations of the Lifshitz formula is to consider the in- 
teraction between two graphene sheets [7j having in mind applications to carbon 
nanotubes. The graphene sheets are described by the two dimensional plasma shell 
model used in [H [9]. For thin shells this model is more appropriate because the 
Casimir or van der Waals forces vanish if the thickness of a dielectric layer goes to 
zero. In [7] the interaction between two such sheets was studied and in [TD] this was 
generalized to the interaction of a material body with a flat graphene sheet and, 
further, to a cylindrical graphene sheet in front of a plane in PFA. 

In the present paper we derive the generalized Lifshitz formula in the first ap- 
proximation beyond PFA for a cylindrical graphene sheet in front of a flat graphene 
sheet or a plane material body with permittivity e(u). In both cases we consider 
two scalar problems corresponding to the TE and TM modes. Regrettably, in this 
configuration the polarizations of the electromagnetic field do not separate (in op- 
posite to the case of a conducting wave guide). However for small separations the 
TM mode dominates and for large separations the ideal conductor case is recovered. 

The paper is organized as follows. In the next section we collect the basic for- 
mulas of the method of [3J and generalize them to semitransparent boundaries. In 
section 3 we calculate the interaction energies in first approximation beyond PFA. 
After the conclusions some appendixes follow with details of the calculations. 

2 Functional integration and semitransparent boun- 
daries 

In this section we start from the method of introducing boundary conditions into a 
functional integral by functional delta functions. The method was used in [11] for the 
calculation of radiative corrections to the Casimir effect. Later it was rediscovered 
[T2] . In [S] and [3] it was used to obtain a finite expression for the Casimir interaction 
energy of a sphere or a cylinder with a plane for Dirichlet and conductor bound- 
ary conditions. In [5] similar results were achieved in a multiscattering approach. 
Here we follow the representation given in [3] and generalize it to semitransparent 
boundaries. 

Before entering the formalism of functional integration we formulate the bound- 
ary conditions which are considered in this paper. There are two types, that on the 
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surface of a dielectric body and that on an infinitely thin plasma shell. We mention 
that both are in fact rather matching conditions relating the fields on both sides of 
the surface. These boundaries are semitransparent in the sense that the correspond- 
ing reflection and transmission coefficients take values somewhere in the interval 
from zero to unity. Both turn into Dirichlet resp. Neumann boundary conditions in 
the limiting case of the frequency parameter going to infinity. 

1. Plasma shell 

Here we consider an infinitely thin sheet filled with a charged fluid in an oppo- 
sitely charged immobile neutralizing background. This is the two dimensionale 
limiting case of three dimensional plasma used for example in the theory of 
metals. It is aimed to describe the 7r-electrons of a graphene sheet or a carbon 
nano tube. For details see [9] and papers cited therein. The interaction of such 
a sheet with the electromagnetic field can be reduced to matching conditions 
across the sheet. For a flat sheet the polarizations separate into TE and TM 
modes and the matching conditions read 

= 0, $V - = 2fi $, (TE) m 

$' + - = 0, = -2 (TM) ' ' 



where 
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Q = 2?r (2) 

m 

is a parameter in parallel to the plasma frequency uj p in the second model. It 
depends on the parameters of the plasma, the density n of the electrons, their 
charge e and mass m. The inverse, is a length which is to be compared 
with the geometric sizes of the interacting bodies. For a carbon nano tube we 
note 1/Q = 1.5/im (see [10], Eq.(ll)), which is by two orders of magnitude 
lager that the corresponding value, say l/u p = 0.02 fim for gold, in the second 
model. The matching conditions for the TE polarization is just the same as 
for a delta function potential of strength Q for a scalar field and up to some 
subtleties the conditions for the TM mode correspond to a potential with the 
derivative of a delta function. As shortcut we will use the subscripts <5te and 
Stm- 

In addition to a flat sheet we consider also a cylindrical one. In that case the 
polarizations do not separate (in opposite to a conducting wave guide). We 
define the analogous scalar problems by the same matching conditions ([!]) as 
in the case of a flat boundary on the radial functions, 



/rib .°° 



(3) 



with polar coordinates (r, (p) in the plane perpendicular the the axis of the 
cylinder, x\ = rcostp, X2 = rsimp. The prime denotes then the radial deriva- 
tive. 
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2. Dielectric body 

We consider a dielectric body with permittivity e filling the half space Xi > 0, 

(e(u) forxi>0, 

In order to consider later the short separation limit we need to include fre- 
quency dispersion. We take the plasma model, 

<«>) = 1 " ^ (5) 

where u p is the plasma frequency. Usually, this model is used for metals. 

Since the surface is flat the polarizations of the electromagnetic field separate 
into transverse electric (TE) and transverse magnetic (TM) modes. The well 
known matching conditions for the corresponding scalar functions (amplitudes) 

®(z) are 

= 0, $' + - = 0, (TE) 
e$ + - $_ = 0, $' + - = 0, (TM) 1 ' 

where $± are the limiting values from the right and from the left of the plane 
and the prime denotes the derivative with respect to x±. Here the Fourier 
transform in the translational invariant directions X\\ = (x 2 ,x 3 ), i.e., in the 
directions parallel to the plane, is assumed, 

-J 6**11*11 $ W)fc|| (a;i). (7) 

The amplitudes satisfy the wave equation 

d 2 \ 

e(z, u) uj 2 -k\+ \ $ Wifc|| (an) = 0. (8) 

In the following we use the shortcuts €te and ctm to identify these models. 

The method of implementing boundary or matching conditions in the functional 
integration rests on the restriction of the integration space in a functional integral 

Z(J) = I ' D<f> l[S(H g [<f>]) exp{-S} (9) 
J s 

representing the generating functional Z (or the partition function) of a field theory 
with a field $ and some action S. Initially, the theory is completely arbitrary and 
it may contain interaction, background fields and boundaries. The functional delta 
function provides then the necessary restriction of the integration space to fields 
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fulfilling additional boundary or matching conditions on a surface S. We assume 
that these conditions can be formulated in terms of a linear functional, 



H z [$] = J H(z,x)${x) dx = , (10) 

where H(z,x) is some integral kernel. For example, Dirichlet boundary conditions 
are given by 

H D (z,x) = S(x-f(z)), (11) 

where the function f(z) describes the surface S = {x\ x = f(z)} and z provides a 
parameterization of the surface S. Assuming the surface S is a plane parallel to the 

(x 2, x 3) -plane at x\ — a in an M 3 . For that a parameterization is z = ( X2 ^ and 

( a \ 

f(z) = I x 2 .In fact, below we need the corresponding quantities for a cylinder of 
\x 3 J 

(Rcos(p 
Rsm ip 

The function H(z, x) for the upper line of the matching conditions (CQ) reads 

H 5 ™(z,x) = 5(x- f(z)) (d n+ - d n _ - 2ti) , (12) 

where d n+ is the normal derivative on the right side of the surface (xi > a) and 
<9 n _ is that on the left side. The function H(z, x) for the lower line of the matching 
conditions ([T]) reads 

90 

H 5 ™ (z, x) = 6(x- f{z)+) -S(x- f{z)J) + —5{x- f{z)) d n . (13) 

us 1 

Here f(z) + assumes to take the limiting value if x approaches the surface from the 
right side. The normal derivative can be taken on either side. Using these functions 
in Eq. ffTUj) one obtains just the matching conditions ([T|) for a plane or, with the 
appropriate choice of z and f(z), for a cylinder. 

It should be mentioned that this method is completely general and that it holds 
for any boundary or matching conditions which can be written in form of a linear 
functional like (jTOj) . Also, it is not restricted to a scalar field. Originally it was used 
in [H] for the electromagnetic field. 

However, there is also a limitation of this method. It cannot been used for the 
matching conditions on the surface of a dielectric body since there are different 
speeds of light on both sides and the condition cannot be expressed in form of a 
delta function in the functional integral. Whether this limitation can be overcome 
is not known at the moment. There exists an attempt to incorporate a dielectric 
body into the functional integral [13] but it looks too complicated and no nontrivial 
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Figure 1: The configuration of a cylinder in front of the plane S 



example was given so far. Also the recent approach of [H] cannot solve this problem 
since it cannot handle the TM mode. 

With the formulation ffTUl) of the boundary conditions the method follows exactly 
the steps given in section 2 of [3]. One arrives at a representation of the Casimir 
energy in the form 



in the space of functions defined on S and the trace is over such functions. D(x,y) 
is the propagator of the initial theory defined by the action S in ([9]) and the trans- 
position in H T (y, z') means that the derivatives in H act to the left. A difference of 
this formula as compared to Eq.(15) in [3] is that we already passed from time de- 
pendency to frequency dependent quantities by means of the corresponding Fourier 
transform. Also, starting from here we work in the Euclidean version. 

The configuration we are interested in is shown in Fig. [1] and we proceed as 
follows. We consider a free scalar field in the presence of a dielectric halfspace at 
X\ > a or a flat plasma shell at x\ = a given by the action S and the functional 
integral (Q without the functional delta function as the initial theory. Further we 
consider this initial theory to interact with a cylindrical plasma shell S with radius 
R. This interaction is incorporated in (Q as matching conditions ([1]) using the 
functional delta function. In this way we formulated two problems, the interaction 
between a cylindrical plasma sheet S and 

1. a flat plasma sheet (shortcuts (5<5)te and (<W)tm), 

2. a dielectric halfspace (shortcuts (eS)^ and (eS)™)- 




(14) 



where K has an integral kernel, 




(15) 
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In the following, to a large extend the formulas are the same for all considered 
problems. The shortcuts will be used for the quantities specific for the considered 
problem. 

In general, other combinations, for example the interaction of a cylindrical shell 
obeying TE conditions with a plane carrying a TM condition, can be considered too 
but we restrict ourselves here to the cases which appear in the physical applications. 

To continue we have to specify the propagator D(x,x') of the initial theory, 
i.e., for a flat plasma shell and for a dielectric halfspace. These are well known 
expressions but we need them is a somewhat specific representation, namely as a 
difference 

D u (x,x , ) = DS\x-x f )-D u (x,x r ), (16) 

where the subscript uj indicates that the Fourier transform in the time variable was 
done. The first term in the r.h.s., D^(x — x'), is the free space propagator for a 
massless scalar field, 

/rl'^h p -ik{x-x') 
JUL - (i7) 

In the following we need the representations 

D^(x-x') = J^ 2 e -^ir*I|) d^fr - si) (18) 

with 

<i„ (*i - <) = < -^P~ (j = y/" 2 + H ) , (19) 
which emerges from fflTl) after integration over k\ and 

(20) 



DS ) {x -x>) = f^J- £ e -^3(,3-4) + ^-,') 4°) 3im (r,r') 

^ ' 171— — OO 



with 



j(o) / f\ I Im(pr)K m (pr') for r' > r , , , \ 



I m (pr')K m (pr) for r > r' 



which is the representation of the propagator in cylindrical coordinates in terms of 
modified Bessel functions. 

In Eq. (fl6l) the addendum D^x^x') describes the boundary dependence of the 
propagator. For the flat plasma sheet we make use of the translational invariance 
in the directions parallel to the sheet and write down a representation in parallel to 

CUED, 

D«(x,x') = J-^J e -^ll-!|) d^( Xl , x [) (22) 
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and similar for Du{x,x'). Then the difference 

j(0) 



4,,fc||(£i,3a) = dl'k u ( x i ~ x 'i) - du,k H (xi, x[) (23) 

is subject to the same boundary or matching conditions as before. 

For simplicity we start from Dirichlet boundary conditions. In that case the 
representation 



holds which can be checked by applying the boundary conditions, i.e., (12"3"1) with 
inserted vanishes if either x\ = a or = a holds. Note that this representation 
holds for the arguments x\ and x[ on any side of the plane at x\ = a. For both 
arguments to the left of the plane the expression simplifies, 

= hi ^ {2a ~ X1 ~ <] < a ) • (25) 

Inserted together with (|T9j) into (I2"3l) and, further, into (TL8I) . just the well known 
propagator with one reflection on the mirror at x\ — a appears. Now we write down 
a similar expression for the the plasma shell. For the TE matching conditions (first 
line in (Op) we obtain 

JSte I -7(|xi-a|+[xi-a|) tcyc\ 

d^i^i) - 2 7 l + 7 /i2 6 ' W 

which can be checked again by applying the matching condition. For both arguments 
to the left of the plane the expression simplifies, 

dX <) = ^ Y^jn e"^ 2 — i) ( Xl , x[ < a) . (27) 



Similar expressions can be found in [15] where the propagator for a delta function 
potential was considered. 

For the TM matching condition (second line in (JT])) we proceed in the same way 
and obtain 

^(zi^i) = Y "jf^ sign^x - a) sign(x; - a) (28) 

which can also be checked by applying the matching condition (note that we are 
here in the Euclidean version and uo 2 enters with the apposite sign as compared to 
Eq.([ID). We note that d^™\xi,x[) has a jump if x\ passes through a whereas the 
derivative is continuous. To the left of the plane the representation simplifies too, 

^(-i, <) = ± e -7(— x[ < a) . (29) 
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Now we need the corresponding formulas for the dielectric half space. These are 
standard and can be found in many places. For the TE mode it holds 

dX(xi,x[) = ±- ^1 {xux[ < a) (go) 

.ii 27 7 + p 



with p = ^Je{iuj)uj 2 + k 2 which is the momentum perpendicular to the plane in the 
medium after rotation to the imaginary axis. For the TM mode we note 

feO^'i) = =i e -<^-<) {xux[ < G ) . (31) 

' 11 27 e(ftj)7 + p 

Obviously, the formulas (1271) . (|29l . ( 1301 and (13T1) have the same dependence on x\ 
and x'± and they can be joined into 

L H {x 1 ,x' l )=d u] , 1 e -T( 2 «-«-*i) (n,^,,) (32) 

with c/^ i7 to be substituted by one of 



J<5te — _L 1 J<5tm - J_ -7^ 

W)7 27 l+7/fi ' w,7 27 o; 2 +7n ' 

Je T E = _L 7jl2 J^TM = nL e(i^)7-p 

w,7 27 7+p ' w,7 27 e(iu))~j+p 



(33) 



Up to a common factor these are just the reflection coefficients of the related scat- 
tering problems. 

Now with the explizite expressions for the propagator D(x,x') of the initial 
theory at hand we return to Eq. (IT4"|) . The next step is to calculate the trace over 
K(z,z'). We chose a convenient basis for that, namely 

\k 3 ,m) = — e ik3XS+im(p (34) 

with respect to which all functions defined on the cylinder S can be expanded. For 
details see [61 [3] . Then Eq. (fT4"j) can be written in the form 

1 f du f dk 3 

E Cas = 2 / 2vr / 2tT trmlnKm > m ' ( 35 ) 

with 

Km tm > = (k 3 ,m I K(z,z') I k 3 ,m') . (36) 

We mention that K m ^ m i is an infinite dimensional matrix labeled by m and ml and 
that In K m>m i is another matrix and that the trace is over the latter matrix. In 
Eq. (j35p the translational invariance along the axis of the cylinder was taken into 
account and -Enc a . fl35|) . is in fact the energy density per unit length of the cylinder. 
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We proceed with dividing K m ^ m i into two parts according to the subdivision in 
Eq.flED and (USD, 

■Km,m' = -K mrn t K mm i (37) 

with 

K t]m> = (h,m\K^\z,z')\h,m'), 

Km,m> = (k 3 ,m \ K(z,z') \ k 3 ,m'). (38) 

Now we remark that the functions H (z, x) appearing in the formulation ffTOl) of the 
matching conditions commute with the averaging in the basis ( l34"j) and we obtain 



K m ,m> = J dr dr' H(r) d W]k3im , m >{r,r') H T (r') (39) 

with 

du,k 3 ,m,m'( r i r ') = (h,m I D u (x,a/) \ k 3 ,m'). (40) 

With the same subdivision as in Eq. ([16]) into free space part and addendum depend- 
ing on the boundary conditions on the plane in x\ = a we define 

4> ) ft 3 ,m ) m'( r > r> ) = rf i°i 3 ,m( r ' r ')$rn,m> ~ diO,k 3 ,m,m'{r, r'), (41) 

where m (r, r') is the free space part introduced in Eq. (T20|) . Obviously its con- 
tribution is diagonal in m and m! . 

Next we need the functions H(z,x), Eqs. (lT2l and ( {TBI for the cylindrical surface 
S. They read 



H s ™( r ) = 5{r - R) (d r+ - d r _ - 2Q) 



2Q 



H d ™(r) = 5(r-(R + 0))-5(r-(R-0)) -5(r-(R + 0)). (42) 

u> 2 



(0) 

K$$P = ^K^™) (43) 



Now we calculate i^ m /- For the TE case we have 



with 



^° )(TE) = j dr dr' H**(r) <i 3 , m (r,r') H s ™ T (r>). (44) 

First we remark that in applying if 5 ™ 7 \r') we have to keep r ^ R so that d® ks m (r, r') 
is continuous at r = R including its derivatives and we come to 

<° ){TE) = Jdr H***{r) $ Mjm {r,R) (-20). (45) 
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Now we apply H TE (r) and we have to take care of the jump of the derivative with 

i (0) ( 



respect to r of d^\ (r, R), see Eq.(J2H). Using I' m (z)K m (z) - I m (z)K' m (z) = 1/z 



we get 

on 

^i° )(TE) = - p- (1 + 2fli2 J m (pi2) ^ m (pi?)) • (46) 
it 

Proceeding in the same way with the TM mode first we define 

<S M) = W^° )(TM) (47) 

with 

k (o)(tm) = J dr dr' H s ™(r) cg^Jr, r') H s ™ T ( r >). ( 48 ) 
Application of H s ™ (r') gives 

K (0)(TM) = F , TM(r) d (0) (ry) £ Ml (4Q) 



-U 2 



and finally we obtain 



^(oxtm) = (l - ^ • (50) 

In this way we calculated the projection of the free space part of the propagator onto 
the cylinder with the matching conditions of the plamsa sheet. These expressions 
alone, i.e., Km ^ or Km U \ inserted into Eq. fl35|) would give the Casimir energy 



of a cylindrical plasma shell (more extactly, since the polarizations of the electro- 
magnetic field do not separate, the Casimir energies of the two scalar problems 
defined by the matching conditions ([T])). This Casimir energy contains ultraviolet 
divergences which, b.t.w., are not investigated yet. 

Now we calculate the second part in (1371) which contains the information on the 
plane in x\ = a. Here, of course, we do not have the cylindrical symmetry. For the 
propagator D(x,x') in K m ,m'i Eqs. (1551) . (fTol) . we use the representations (I2"2"j) with 
(I32"|) . We note that D(x,x') and its derivatives are continuous at the cylinder (the 
jumps are on the plane at X\ = a). Therefore, in each function H(r), Eq. (1421) . only 
the last term contributes. In this way we come to 

H 5te D(x,x')H 5teT = (-2fi) 2 D{x,x')\ (51) 

\r=r'=R 

and 

— J d r d r ,b(x,x')\ (52) 

IjJ 1 J \ r =r'=R 

with the already mentioned cylindrical coordinates x± = rcosip, X2 = rsimp, x[ = 
r' cos if', x' 2 = r'sinip'. 
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We have to take these expression in the basis ([Mil and obtain with ([22]) and (132|) 
#S = fem|(-2fl) 2 5(^'), |fc 3 ,m'> 

lr = r' = fl 

= d£; 2 / m /*, e" 2 ^4 7 (53) 



with 

/L = / 

2tt 







Introducing new variables 



&2 = P sinh 0, 

7 = pcosh^, (55) 



the transforming the integration over ip into an integral representation of the mod- 
ified Bessel function, the f m become 

fm = e me I m (pr)- (56) 

With this the final representation of K m m , is 

#£3 = 4fi2 l m(pR) I m ,(pR) K, m+m , (57) 

with 

/■oo 

JC m+m/ = d6 cosh ((m + m')6) 2p cosh 9 d u „ e - 2apcoshe . (58) 
</o 

We note that in the limiting case of Q — > oo where the matching conditions ([I]) 
turn into Dirichlet boundary conditions, 2pcosh6 l d Wi7 — > 1 holds and /C m+m / ( 158|) 
becomes an integral representation of the modified Bessel function, 

fcm+m' — * K m+m i(2ap), (59) 

which coincides with the corresponding formulas in |6] and [3j. However, for finite 
Q we are left with the integral representation ( 1581 ). 

For the TM case we proceed in the same way. After applying the functions H(r) 
( 152]) and projecting on the basis | k 3 ,m) we get 

= d r d r ,(k 3 ,m I D(x,x') l I k 3 ,m') . (60) 

m,m \UJ Z J \r=r'=R 

The matrix elements can be calculated in parallel to the TE case. The radial deriva- 
tives apply only to the Bessel functions I m (pr) since /C m+m / does not carry any 
dependence on r and we end up with 

k™, = UpR)C(pR) £m +m '. (61) 
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Again, in the limiting case -> oo we reobtain the corresponding expression for 
hard boundary conditions, Neumann ones in this case. 

By means of formulas ( l46l) , (|50l) and (|57l) , ( J6TT) we calculated the matrixes entering 
K m ^ m i in Eq. fl37|) . Since K^ m , is diagonal we can easily separate it and the logarithm 
in Eq.(j2SD becomes 

lnK m>m , = In (K$5 m , m ,) + In (S mjm > - A m ^ m ,) (62) 

with 

In (1621) the first term gives the Casimir energy of the plasma shell cylinder alone. The 
decisive point of the method is now that this part carries all ultraviolet divergences 
which are in -Ecas, Eq. lfUL and that it does not depend on the distance between 
the plane and the cylinder. Therefore the distance dependence is solely contained 
in the second term in the r.h.s. of (1621) . It gives rise to the distance dependent part 
of the Casimir energy, 



1 f duj f dk 3 

-C^Cas = 77 / 7; — / 7; — t r m 111 {Om,,m' ~ A. m r 

2 J 271 J ATT 



(64) 



which does not contain any ultraviolet divergence. It serves as the basic represen- 
tation for the further elaboration. 

Inserting (14"6"|) . (I50|) and (157j) . (l5"Tp we simplify the expressions to some extend, 

< E m '=^ TE £m+m> (65) 

with 

^ TE = 1 (66) 

1 + 2<nRI m {pK)K m (pR) K ™{pR) 



and 
with 



<™ = £m +m > (67) 
D<5tm _ \ I'm'(P R ) (ao\ 

~ 1 - 1 K> (pRY { } 

1 2UR P 2 I' m (pR)K' m (pR) ™\P V 

In this representation the information on the cylindrical plasma shell is in the first 
factors on the r.h.s., i.e., in R Ste and in R s ™, whereas the K m+m i by means of 
Eqs.f l58p and ( 1321) carry the information on the plane at x\ = a, whose different 
cases will be indicated in the following by a corresponding superscript, JC^f m , or 
fcm+m> f° r example. 
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3 Small separation and first correction beyond 
PFA 



In this section we generalize the small separation expansion of section IV in [3] to the 
semitransparent boundaries of a plasma sheet and a dielectric half space interacting 
with a cylindrical plasma sheet. For the dielectric halfspace with dispersion e(a>) we 
restrict ourselves now to the plasma model, i.e., to a permittivity given by Eq.Q. 
We start from the representation ( }64"1) of the interaction Casimir energy with the 
matrix elements A^, and A 5 ^^, given by Eqs. ( |65j) and ( 1671) . The idea of the small 
separation expansion is is that only high momenta, i.e., large all, u, k%, m and m', 
contribute and that one can use the asymptotic expansions of the matrix elements. 
Further, it is useful to expand the logarithm although at the end it needs to be 
summed up again. We note that this asymptotic expansion is the same as done for 
hard boundary conditions in [3], however the expressions depend in addition on the 
parameters Q and uj p and they are more involved. 
We start by expanding the logarithm in Eq.flBi 



-1 f°° clu f°° dk 3 1 f°° f°° f°° 

with 

The next step is to make the rescaling u — > uo/R and k% — > k 3 /R which makes R 
disappear everywhere except for a factor 1/R 2 in front of the whole expression and 
for the exponential in A m+m /, Eq. fl58|) . where it combines with a into 



l = ^ = 1 + £ - < 71 > 

where we introduced the distance L between the cylinder and the plane and e = j|, 
which is the small parameter at short separation. The other two length parameters in 
the considered problem are 1/Q and l/u p . Of interest are values in the nanometer 
and micrometer region. The short separation expansion implies 1/Q << R and 
l/ui p « R. However, with respect to the separation L we assume them to be 
of the same order, i.,e., 1/Q ~ L and l/u p ~ L. As a consequence we have to 
consider the combinations QRe = QL = Ql and u) p Re = uj p L = ujl which will 
appear in the following as being of order of one. The latter notations will be used 
in the intermediate steps. In order to shorten the notations, in the following we 
will put R — 1. The correct dimension of the energy can be restored by -Ecas ~~ > 
E C as/R 2 . In the final expressions, only the parameters VLl an d wl will appear and 
their dimensions are restored by fi^ — ► QL and lvl ~~ > uj p L. 

We proceed with the following substitutions of variables in f lBT?]) . Taking into 
account that the integrand is an even function of u, ks and m we first substitute 
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io = p sin a and = p cos a and then 



p = -y/l — r 2 , m = -r, rii — > rii\J — (72) 

after which the Casimir energy can be written in the form 

-e~ 3 ^ 1 f°°dt , f w/2 da f 1 , 
27r s + 1 Jo t Jo n/2 J 

m ... n s - .M , (73) 



oo •/ — oo 



where .M is still given by Eq. (170l) . We note that in the TE case Ai does not depend 
on a and r and these integrations give a factor of unity. However, in the TM case 
M. depends on both these variables. 

In Eq. (1731) we already made the first step in the asymptotic expansion by substi- 
tuting the summations by integrations. The next step is to expand in A4 the A m m /, 
which are given by Eqs. (165|) and ( 1671) . We start from the factors R 5te and R s ™. 
Here we have simply to insert the uniform asymptotic expansions of the modified 
Bessel functions. We note that the exponential factors are the same as for hard 
boundaries since they cancel in the (^-dependent denominators. We obtain 



i^TE „ r (5 TE ) ( X + p(«TE) ^ + Q(5 TE ) e + \ t 74 \ 

e^O IT " 

and 



rStm r (5 TU ) ( 1 + p(5 TM ) v ^ + g(5TM) e+ _ ) (75) 

with 



e^O 7T 



i + tAV 



1 + ty l /U L 



where we introduced the notation 



Vl — t 2 sin a. (77) 



Details on the calculation and explicit expressions for P\ St *> and Q^ St *' are given 
in the Appendix A, Eqs. (IA.3l) . The factor r/o in the exponential will be considered 
below together with the corresponding one from JC m+m f. 

We continue with the expansion of JC m+m i. In the case of hard boundary condi- 
tions this was due to (|59p simply the asymptotic expansion of the Bessel function. 
Here we have to consider the integral representation (1581) . However, since it is very 
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close to a known integral representation of the Bessel function we can use standard 
methods, namely a saddle point expansion. This done in Appendix B where also the 
result is stated. It consists of a factor 0, Eqs.( 1B.9l )- (1B.lll) . which is different for the 
models considered and factors which are the same as for hard boundary conditions 
in [3]. These factors are collected in ip, Eq. (1B.13l) . With all these factors at hand 
we obtain for the functions A m+n , m+n i entering Ai, flTUl) . 

A m+ n,m + n' = ^ r^r^ (l + ^ a^ff + e a% + ...), (78) 

where the factor in the exponential is the same as in [3J, Eq.(57), 

rf* = 2t+ {n-n'f. (79) 

Formula flTHj) has the same structure as Eq.(56) in [3] with the only differences that 
the functions a^ffl and a^ n , are different (they are displayed in Appendix C) and 
that the additional factors r^ cy1 ^ and r^ plane ^ are present. These are in the given 
variables just the reflection coefficients for the corresponding boundary conditions 
on planes whereby r^ 1 ' originates from R 6 ™ in (|65l) and R s ™ in ( IBTl) . i.e., from 
the cylinder, and r^ plane ^ originates from in A m+m >, (1581) . i.e., from the plane. 
There are four expressions, 

for plasma sheet, TE mode, 
for plasma sheet, TM mode, 



for dielectric, TE mode, 



for dielectric, TM mode, 

whereby r^ cy1 ^ can be equal to the first two and r( plane ' ) to any of them. 

In (1811 we used the notation y, Eq. (l77|) . The variables r and a enter the reflection 
coefficients just in this combination. In the leading order for small e this is the only 
dependence on r and a and one of these integrations can be carried out even in the 
TM case. In higher orders in e, from the expressions a^) and a^ n ,, a dependence 
on both variables comes in. However, the dependence on r is polynomial and can 
be cared for by means of the simple formula 

r^To [ 1 drr^f(VT^sina)= f ' dy (1 - y*Y f(v) ■ (81) 
Jo n/2J ^7i\.(n + l)J 

Next we insert ( 1781) into A4, ( I7U1) . and after a reexpansion we obtain 

3+1 

J\4 = ( 6 ^ 2 g-2(s+l)t-m r (cyl) r (plane) ^as 



,(cyl, plane) 



ty 2 +Q L 



coi +t 



(t+y/L 
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with 



M" = l + v^E<S +1 (82) 



i=0 



+£ J ^2 a n/,J i+1 a nj(n j+1 + a nln i+1 ) + • • • 
\0<i<j<s i=0 



and 



Vi = ^2 (ni - n i+1 ) 2 , (83) 

i=0 

where, in order to include all contributions into the sum signs, we have to put 
n = n s+ i = formally. With these formulas the Casimir energy Ecas, ( |73l) . becomes 

° as 2n ^s + lj t yfa Jo k/2 L 



=o 







^drn f™^ ^ ^ (cyl ) r(plane )^+i ^as (g4) 



7T 

This expression for the Casimir energy is to some extend the final formula for the 
short separation expansion because the remaining integrations cannot be carried 
out in an explicit form (except for that over the ni). It represents the generalized 
Lifshitz formula beyond PEA. In the following we discuss its basic features. 

First of all, the dimensions are restored by multiplying the whole expression 

be R~ 2 . With e~ 5 ^ 2 R~ 2 = j?\pt the dimensional factor known from the PFA is 
restored. The dimensions inside the above expressions are restored by means of 
Ql — > VtL and Ul — > uj p L. Further, it is obvious that the limes of hard boundary 
conditions is recovered for Q — > oo and uo p — > oo since we arrive just to the same 
expressions as in [3]. 

The first nontrivial result which follows from is a confirmation of the PFA 
for semitransparent boundaries derived in [10]. It comes about from the leading 
order in e, i.e., with Ai as — > 1 in ( 1841) . With this, there is no dependence on the 
besides in r]i and the integrations over the rii can be carried out using Eq.(66) in 
[3] delivering a factor of (s + 1) _1//2 . Next we observe that there is no dependence 
on r so that formula (IHTj) can be used with n = 0. Finally, the summation over s is 
expressed in terms of a polylogarithm, Li s (z) = Y2^=i z n /n s , and we arrive at 



E Z A = ^fik? V * [ * ' 3/2 / * [WW «"*] (85) 



for one of the scalar problems with r^ cyl ) and r^ plane ^ given by Eqs. flHOj) . Up to 
differences in notations this coincides with Eq.(33) in [TO] . 
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In order to represent the result in a more instructive manner we rewrite (f8"4"|) in 
the form 



E, 



7T" 



Cas 



(f (QL, o; p L) + | A h Wp L) + . . .V (86) 



1920^^2 V L 

where the factor in front is the Casimir energy for hard boundary conditions (Dirich- 
let for TE and and Neumann for TM, both give the same) in PFA. The function 
/o represents the relative decrease of the Casimir energy in the considered model in 
PFA. The function fx is the first contribution beyond PFA. It is also written relative 
to the hard boundary case for which we have [3] 

for TE, 

(87) 



36 



for TM. 



Let us consider the PFA, i.e., the function / , in more detail. It follows from Eq.(l85 
with r^ 1 ) and r <plane * ) for the models formulated in section 2, 



/ ((M)TE (ftL) 


480^ 

7T 9 / 2 


# s) ™(nL) 


480^ 

7T 9 / 2 


/f )TE (fiL,^L) 


480^ 

7T 9 /2 


ft ] ™m,u; p L) 


480^ 

7T 9 /2 



dt t 3/2 Li 
dt t 3 ' 2 
dt t 3/2 Li 



3/2 
1 



exp(-2t) 
_{l + t/QL) 2 _ 

exp(-2i) 



dy Li 



3/2 



dt t 3/2 / dy 



3/2 
1 



l + fy 2 /OL) 2 
^(ujpL) 2 + t 2 -t exp(-2t) 
t + y/(u p L) 2 +¥ 1+t/QL 







Li3/2 



t (t - v/KL)2+t2j y 2 + ( ^ L) 2 exp( _ 2t) 



,L) 2 + t(t+ y/(iu p L) 2 + i 2 ) y 2 l + ty 2 /QL 



oo, o) p L 



1 for all four functions. In that case L<<^,L<< — and L « R 



The limiting value for hard boundary conditions which appears for VtL 
oo is f 
holds. 

The opposite limiting case is^<<L<<i? and -j- « L « R. It corre- 
sponds to small separation as compared with the plasma wave length. This is the 
nonretarded regime or the case when the Casimir forces turn into the van der Waals 
forces. It was already mentioned in |10j . 

To obtain that limiting case, for the TE cases one simply expands the polyloga- 
rithm formally for small L and the remaining integration then gives 



fo 



/ Wte (OL) 
{e5)TE (nL,uj p L) 



QL-*0 



s 4 mr . 



7T J 



420 2 
— r v p \lL 

7T* 



59) 
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Figure 2: The function f( eS ) (x) appearing in ( l90l ) 



In this limit the TE mode is subleading. The leading contribution comes from the 
TM mode. Here it is impossible to simply expand the arguments of the polyloga- 
rithm because then the remaining integrations become singular. Instead on has to 
make the substitution either y — > yuj p L or y — > y\/VLL after which the limit L — > 
can be taken, 



/| W) ™(fiL) 



840^ 



7T 



9/2 



(5<5) 



/ (e5) ™(fiL,^L) 



840^ 

r 9/2 



W>(f) 



(90) 



with 



/ 



(55) 



and 



./ 



(e<5) 



dt t 3/2 / Li 3/2 
'o Jo 
0.254, 



/* oo roo 

/ t 3 / 2 / dy Li 3/2 
Jo Jo 



exp(-2t) 
(i + ty 2 ) 2 



exp(-2t) 



1.39. 



l+t?/ 2 W 2 L/fi)(l+2tVy 

is shown in FigfSJ Its limiting values are (x) ~ 



z^O 



(91) 



(92) 



1.62Vx 



The function /( e<5 ) 
and (x) 

x— >oo 

The functions /o for the (M)-case, i.e., for the interaction of the cylindrical 
plasma sheet with a flat plasma sheet, are shown in FigfSJ For large argument 
they go to unity which is the hard boundary limit. For smaller argument when 
the sheets become more transparent they decrease whereby the function for the TE 
case decreases faster. The behavior at the origin is given by the upper lines in 
Eqs.f l89p and (!90|) . The corresponding functions for the (e<5)-model depend on two 
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Figure 3: The function fQ S (QL) representing the decrease of the Casimir energy 
at small separation relative to the hard boundary case for the TE and for the TM 
modes in the interaction of the plasma shell cylinder with a flat plasma shell in PFA. 




Figure 4: The function f^ S) ™(nL, oo p L) with u p = y/QX/L for A = 1(F from j = -2 
(lower curve) till j = 2 (upper curve) in equal steps. 



parameters. Their general behavior is similar to the (55)-case and as an example 
we show in FigJUthe function fQ tS ^™(QL, u p L) for several values of the ratio of its 
arguments. Also for these function the limiting value for large arguments, i.e., for 
hard boundaries is unity and the limit for small arguments is given by the lower 
lines in Eqs.flHH]) and flip. 

The first contribution beyond PFA is given by the functions fi in Eq. (|86p . These 
appear from the contributions proportional to yfe and to e in the representation (|84|) 
of the energy with ajffi and a® n „ Eqs.flCTIl till fl08|) . inserted into M as , ([82]). The 
next step is then to carry out the integrations over the rij. These integrations are 
Gaussian since the exponential rji is quadratic in the rij and Ai as is polynomial. 
As seen from the explicit formulas, is odd in the rij hence the contribution 

proportional to \ft vanishes. So we are left with the integrations of the terms 
proportional to e. These integrations can be carried out in the same manner as in 
[3|. However, the results are too big as to be displayed here. After that formula 
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Figure 5: The function /q TB (OL) and f[ SS ^ TE (QL) (lower two curves) and their 
ratio f( ss > TE (QL) / f^ TE (QL) (upper curve) for the TE-mode in the (<W)-model. 

([51]) can be used since the appearing expressions are polynomial in r. Again, the 
integration over y remaining in (1HTT) can be carried out trivially in the TE case but 
not in the TM case. The next step is to rewrite the summation over s in terms 
of polylogarithms. After that the expressions for the functions f\ become simpler. 
They are represented by formulas similar to that for the fo, Eqs. (1881) . however more 
lengthy and for the (<5<5)-model the explicit expressions are displayed in Appendix 
D, Eq.(EI]). 



It is interesting to note that the functions fo and f\ are numerically quite close 
to each other. For the TE-mode in the (<W)-model both functions are shown in 
Fig.©. Only for small values of the argument QL there is a significant difference. 
It comes about because in this case it is impossible to expand the integrand for the 
function for fi simply in powers of QL. Instead one has to make a rescaling of the 
integration variable and the behavior comes out to be 

/W-dlL)^™^ mf\ (93) 

This behavior is different from the corresponding one in PFA (upper line in Eq. (j89l) ). 

For the TM-mode in the (M)-model both functions are shown in Fig. (J6j). Here 
the behavior for small argument is the same as in PFA and can be calculated by the 
same substitution. The result is 

r(H)i 



~ 0.92 VOL. (94) 

As a consequence, and as seen from the upper curve in Fig. (J6J), the ratio \ S s) TE — 

of the two curves takes finite values for all QL. We mention that like in PFA the 
TM mode dominates at small separation. 

Finally we discuss the results for the (e<5)-model. Again, the functions f\ are 
very close to the functions fo m PFA. So we restrict us to show in FigJ7]the function 
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Figure 6: The function / (QL) and /} (QL) (lower two curves) and their 
ratio f^ TE (^lL)/flf S ^ TE (^lL) (upper curve) for the TM-mode in the (<5<5)-model. 
The function f^ TE (QL) takes always the larger values as compared to f^ TE (fiL). 




QL 



Figure 7: The function / 1 (e5) ™(fiL, uj p L) with u p = y/QX/L for A from A = 0.01 
(lower curve) till A = 100 (upper curve). 
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5 10 15 20 

Fi gure 8 : The ratio / 1 (£<5) ™(fiL,c; p L)// (e,5) ™(fiL,c ( j p L) in the (e5)-model with oo p = 
y^HX/L for A from A = 0.01 (lower curve) till A = 100 (upper curve). 

j(t6)TM ujpL) for several values of the ratio of their arguments. It is more in- 
structive to consider the ratio of the functions, TM (OL, u p L) / fff^™ (Q.L, oj p L), 
in FigJSJ Again, the function fi is always a bit larger than the function f keeping 
a finite ratio. 

4 Conclusions 

In the forgoing sections we have derived the generalization of the Lifshitz formula for 
the interaction of a cylindrical plasma sheet with a flat plasma sheet ((<W)-model) 
and with a dielectric half space ((e<5)-model) in PFA and in the first approximation 
beyond. The limiting values for hard boundary conditions and for a separation small 
as compared with the plasma wave lengths of the model (which corresponds to the 
nonretarded regime) are considered in detail. The previously known results were 
reproduced. The results in PFA and beyond are represented in form of functions fo 
and fi describing the deviation from the hard boundary limes, numerically computed 
and represented in a number of figures. In general words the outcome is that these 
functions are very close to each other, i.e., the functions describing the ratios of the 
two functions fi and fo are close to unity. This allows a simpler representation of 
the result for the needs of practical calculations, namely as a pure superposition of 
PFA and the first correction beyond (as calculated so far in [5] and [I]) and the 
reduction function taking care of the frequency dispersion (or transparency) of the 
boundaries which can be calculated for flat boundaries. In terms of formulas that 
would imply in a representation like flHBl simply to take the function fo in place of 
the function f\. It can be expected that this property holds also if one considers for 
example a dielectric medium with a dispersion given by optical data. 

Open questions are still an extension of the method to a sphere in front of a 
plane and for accounting for the non separation of the modes of the electromagnetic 
field. However, for small separations the TM mode will prevail anyway and for large 
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separations which is equivalent to the hard boundary limes both modes are additive 
(for the wave guide geometry). So it can be expected that for practical calculations 
the above results will be sufficient at present time. 
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Appendix A 

In this appendix we calculate the uniform asymptotic expansions of the factors R 5te 
and R s ™ appearing in Eqs. (1741) and ( 1751) . We insert the well know uniform asymp- 
totic expansion of the modified Bessel functions into these formulas. The exponential 
factors in the ^-dependent denominators cancel and that from the quotient of Bessel 
functions will be considered later combined with that from the JC m+m >. For the mo- 
ment we denote then by 770- We introduce the intermediate notations z = p/m and 
z' = p/m! and obtain 



i° 1 [m /l + z /2 V / u^z) u^z') 



R^ = - ^== X -At-^) 1 + ^ + ^1 + ... , (A.1 



7T 1 i/mV \ m' \1 + z 2 J \ m 

where Ui(z) = (3(l + z 2 ) — 5)/(24(l + z 2 ) 3 / 2 ) results from the first Debye polynomial. 
In the same way we obtain for the TM case 



7T H f V m' \ 1 + z' 2 J V z \ m m' 

(A.2) 

with v^z) = (-9(1 + z 2 ) + 7)/(24(l + z 2 f /2 ). Now by make the substitution ([72]) 
and expansion with respect to e delivers the expressions stated in Eq. flT4"|) and fl75|) 
with 

p5TE = (n(t-n L ) + n'(n L + t))r 
y/i(n L +t) 

q5te = m^( { - 6{3r2 - 2)n2 - 12n ' T2n - 5T2 

+6n' 2 (5r 2 - 2) + 3) Q 2 L + It (6 (5r 2 - 2) n' 2 
- (I8n 2 + 5) r 2 + 3) tt L + t 2 (6 (5r 2 - 2) n 2 + 12nV 

-5r 2 + 6n' 2 (5r 2 - 2) + 3) ) , 



/ 2 
- n 
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pSTM r (nty 2 + n'ty 2 - nVl L + n'Vt L ) 

Vt(ty 2 + n L ) 

Q s ™ = * ( t 2 (-6 (3r 2 - 2) n 2 + \2n'T 2 n + 7r 2 

-6n' 2 (3r 2 - 2) - 9) y 4 - 2VL L t (6 (3r 2 - 2) n 2 

+ (I8n 2 - 7) r 2 + 9) y 2 + U? L (6 (5r 2 - 2) n 2 - 12n'r 2 n 

+7r 2 -6n' 2 (3r 2 -2)-9)). (A.3) 
We took into account that Qe = Q L and oj p e =ujl must not be considered as small. 



Appendix B 

In this appendix we calculate the asymptotic expansion of the expressions 1C m+m > in 
(155]1 and (157)) . For this we use a saddle point expansion. Let 

poo 

K= d6 g{9) e~ Xh{6) (B.l) 



be the integral to be considered and A a big parameter. Assuming the function h(6) 
has a minimum for some 8 e (0, oo), the expansion for A — > oo is 



/ 5h 2 3 h 4 \ hg 1 g 2 \ 1 

I ;r — — On — — + 



with 

* = g ° + [{^-^J g °-M + ^)x + ---' (B ' 3) 

where gi = g^\9o) and hi = h^(9o) are the derivatives in the minimum, h'(9o) = 
and r](z) is known from the uniform asymptotic expansion of the modified Bessel 
functions. 

For applying this expansion to K. m+m >, we have first to rewrite Eq. (1581) in the 
form 

JC m+m , = - d9 2pCOSh# d u „ e -2apcosh0-(m+m')fl ( R4 ) 

2 Jo 



and then to identify 



g{6) = 2pcosh6 l e? w/y , 

/i(0) = e (2ap cosh + (m + m')0) . (B.5) 



With the intermediate notations 



\ = m + m > z = 2p{l ± e ) (B.6) 
m + m 
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the expansion reads 



7T 




-\ V (z) 



2A (l + ^ 2 ) 1 /4 



+ 



A 



5-3(l + z 2 



9o 



:9o + 



(B.7) 



9i 



+ 



#2 



24(l + z 2 ) 3 / 2yu 2(1 + z 2 ) 2(1 + Z 2 ) 1 / 



+ 



We remark that the function h(0) has just the same structure as in the known 
integral representation of the modified Bessel function and that the whole difference 
is in the function g{6). 

Next we insert the substitution (j72p and reexpand taking into account that 
K, m+m i enters M. in flTUl) in fact with m —>■ m + n and m! —>■ m + n' . Using 
( 132]) for in we obtain 



m+m 1 




Tie 1 
it T 



(B.8) 



where for in dependence on the case considered one has to insert one of the 
following, 



(n + n')n L y/ir 



n L + t (n L + t) 2 
i 



(B.9) 



(p L ((5r 2 - 3) ttl + 2t (12 (r 2 - l) n 2 + 24n' (r 2 - l) 



n 



n 



48t(fi L + 1) 3 

+24tr 2 + llr 2 + I2n' 2 (r 2 - l) - 9) tt L + t 2 (24 (3r 2 - l) 
+48n' (3r 2 - 1) n + 48tr 2 + 41r 2 + 24n' 2 (3r 2 - l) - 15)) )e + O (e 3/2 ) 
VL L (n + n')n L VtTy 2 



ty 2 + n L 
1 



(ty 2 + n L y 

3 (Q L (t 2 (24 (r 2 - l) n 2 + 48n' (r 2 - l) n + 48tr 2 + 7r 2 



(B.10) 



48t (ty 2 + fi L 
+24n' 2 (r 2 - 1) - 9) y 4 + 2Vl L t (12 (3r 2 - l) 
+24n' (3r 2 - l) n + 24tr 2 + 13r 2 + 12n' 2 (3r 2 - l) - 3) y 

+Ql (3-5r 2 ))) e + 0(e 3 / 2 ), 

2 ((n + n')ujlVtr) 



Vjjjf + t 2 - t 

' (cu| + t (-tr 2 - 2 v ^|Tt 2 r 2 + t 

2K + t 2 ) 3/2 (t+ y^TT^y 



(B.ll) 
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((((- 2 -iK 



+t (JujI + V (3r 2 - 1) + 1 (5r 2 - 2)) cu 2 + i 3 + y/u* + t 2 ) (4r 2 - l)) 
+2n' ((r 2 - 1) ui + t (jfi4 + 1? (3r 2 - l) + t (5r 2 -2)^1 



n 



+t J | t + \lu)\ 



t 2 ^j (4r 2 - n + 2t (uj 2 l + t 2 ) (lo 2 l + t (t 



ool + t 2 ) r 



0' 



+n' 2 ((r 2 - 1) + t (yfi% + t 2 (3r 2 - l) + t (5r 2 - 2)) u 2 L 

+t 3 (t + yf^j (4r 2 -l)))cl) 

(v^IT^-t)r 2 (v^iT^-t)(5r 2 -3)\ 

-| ^ ' y ± ' e + O (e 3/2 ) 

2(u 2 +t(t+ v^i + t 2 )) 48t(t+ v^I+t 2 ) / 

*(>/"2 + * 2 -*) ^"^ 2((n + ^V 2 y^ry 2 K + tV)) 

7 \ 2V e l"--'-^/ 



X 



48 (w£ + t (t + v^J+f 5 ) y 2 

24c 2 r 2 K + t 2 y 2 ) (u 2 L +t(t+ ^JjT^) y 2 ) y' 

' fa K + t¥) ((1 - 2r 2 ) W i + t (yQ^V 



K + 1 2 ) 



+t (- (2y 2 + 3) r 2 + y 2 + 1)) u, 2 - (t + y^ 2 + t 2 ) (3r 2 - l) y 2 ) y 2 ) 

+ j-r^f ( 48 -i K + *V) ((((V + 3v^i+^) - 2 - ty 2 



u, 2 +tz)u,i + t 2 i(5ty 2 + Ju, 2 



t 2 (3y 2 + 4)^ r 2 -2ty< 

-y/<4 + * (y 2 + 1)) "2 + (t + \M + * 2 ) ( 4 ^ 2 - 1) ^ 2 ) ' 

+2n' (((V + SyffrH*) r 2 - ty 2 - y^ 2 + u> A L 
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+t 2 ((sty 2 + + t 2 (3y 2 + 4)) r 2 - 2t</ 2 - + i 2 " (y 2 + l)) 
+t 4 (t + v^I + * 2 ) (4r 2 - 1) y 2 ) n + 2t K + **) ^ ((*V 



2 + Ju* + i 2 u, 2 



L 



+t 2 ( t + ^ 2 + t 2 j y 2 j + n' 2 f f f tj? + Sy^ 2 + t 2 \ r 2 - ty 2 
^ ui + t 2 ((sty 2 + y/ul + t 2 (3y 2 + 4)) r 2 - 2ty 2 
+ ? + 1)) <4 +t 4 (t+ sp L + t 2 ) (4r 2 - 1) /)) y^ 

-^f^ (<4 + * (* - \fti + *) y 2 ) (f L + t(t+ \pi^ 2 ) v 2 ) 

+0 (e 3 / 2 ) . 

These expansions were done machined, of course. 

Finally we collect the contributions resulting from the prefactor in (IB. 21) together 
with the expansion of the exponentials in JC m+m ' and in R Ste and R s ™ (which was 
denoted by r] in (lA.ip and (IA.2D ). Up to a common factor these collect into 

(n + n') (2n 2 - An'n + 2n' 2 - At - 1) r n 

* s 1 + 1 " ^ 

+— (-12(n + n') 2 (n 2 - 2n'n + n' 2 - 2t) r 2 + 3 ((5r 2 - 2) n 2 



n 4 



+2n' (5r 2 - 2) n + At (r 2 - l) + n' 2 (5r 2 - 2)) + 2 (-7 (3r 2 - l) 
+An' (3r 2 - l) n 3 + 6 ((3r 2 - l) n' 2 + 2t (r 2 - l)) n 2 
+An' ((3r 2 - l) n' 2 + 6t (r 2 - lj) n 

+6(n + n') 2 (n 2 - 2n'n + n' 2 - 2t) 2 r 2 + 12t 2 r 2 + 12n' 2 t (r 2 - l) 
-7n /4 (3r 2 -l)))e + 0(e 3 / 2 ) (B.13) 



and are used in formula (1781. 



Appendix C 

Here we collect the expressions for a^ 1//2 - ) and a^ 1 -* in Eq. (j78j) . These are the gener- 
alizations of Eqs. (B13) in [5] and for Q L — > oo and cu^ — > oo the old expressions 
reappear. The explicit formulas are quite lengthy, 



= — ^- ( - (-2(fi L + t)n 3 + 2n'(fi L + t)n 2 + (fi £ (2n' 2 + 4t - l) 



(1/2) = r 
(<w)te " 2v^(fi L + t) 



+t (2n' 2 + At + 5)) n + n' (0 L (-2n' 2 + 4t + 3) + t (~2n' 2 + At + 5))) J, (C.l) 
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— — i ^ f 48r 2 t 3 - 72nV - 72n /2 t 2 + 192nVt 2 + 192n ,2 r 2 t 2 

48t(fi L + 1) 2 V 

+ 192rm'r 2 t 2 + 48Q L r 2 t 2 + 21rV - 48nn't 2 - 3t 2 + 2An 2 Q L r 2 t 
+26An' 2 Q L r 2 t - 8 (I8n 2 + 5) Q L r 2 t + A8nn'Q L r 2 t + 22Q L T 2 t 
-2An 2 Q L t - 120n /2 fi L t - 48nn'fi L t + 6Q L t + A8n 2 Q 2 L - A8n' 2 Q 2 L + 9Q 2 L 
-72n 2 Q 2 L r 2 + 120n' 2 Q 2 L r 2 - A8nn'Q 2 L r 2 - 15fi|r 2 

+ (fi L + t) 2 (-84rV + 28n 4 + 48n rV - 16n n 3 - 24n'V + 72n' 2 rV 

+48trV + 30r 2 n 2 - 48tn 2 - 12n 2 - lQn' 3 n + A8n' 3 r 2 n + Q0n'r 2 n 

+96n'tr 2 n - 2An'n - 96n'tn + 28n' 4 - 12n' 2 - 84n' 4 r 2 + 30n ,2 r 2 

+48t 2 r 2 + 48n' 2 tr 2 + 24tr 2 - 48n' 2 t - 2At) 

+(n + n') (-24n(* - fi L ) (H L (2n 2 - 4n'n + 2n' 2 - 4t - l) 

+ (2n 2 - An'n + 2n 2 - At - 3) *) r 2 - 24n'(fi L + i) (fi L (2n 2 - Ann 

+2n' 2 - At - 1) + (2n 2 - 4n'n + 2n' 2 - At - 3) t) r 2 ) (C.2) 

+ (n + n') 2 (p L + tf (2A (n 2 - 2n'n + n' 2 - 2tf r 2 

-24 (n 2 - 2n'n + n' 2 - 2t) r 2 ) - 24 (2n 2 - An'n + 2n' 2 - At - I) t{tt L + t)r 2 ) ) , 

— F — fn £ (2n 3 - 2n'n 2 - 2n' 2 n - Atn - 3n + 2n' 3 + n - Ant) t 

-{n + n')t (-2n 2 + An'n - 2n' 2 + At - 3) ry 2 ^ , (C.3) 
24(n + n') (nty 2 + n'ty 2 - nVl L + n'Q L ) 



n 2 



48 \t(ty 2 + n L ) 2 

((2n 2 - An'n + 2n' 2 - At + l) ty 2 + fi L (2n 2 - 4n'n + 2n' 2 - At - l)) r 2 ) 

+ = 1 - (24(n + n') 2 (2n 2 - 4n'n + 2n' 2 - 4f - l) y 2 r 2 ) 
ty l + iZ £ V ' ) 

+ i (^2 (-12(n + n') 2 (n 2 - 2n'n + n' 2 - 2t) r 2 + 3 ((5r 2 - 2) 

+2n' (5r 2 - 2) n + At (r 2 - l) + n' 2 (5r 2 - 2)) 
+2 (-7 (3r 2 - 1) n 4 + 4n' (3r 2 - l) n 3 + 6 ((3r 2 - l) n' 2 
+2t (r 2 - 1)) n 2 + 4ra' ((3r 2 - l) n' 2 + 6* (r 2 - l)) n 

+6(n + n') 2 (n 2 - 2n'n + n' 2 - 2tf r 2 + \2t 2 r 2 
+I2n' 2 t (r 2 -l) -7n' 4 (3r 2 -l))) 



t (ty 2 + n L 



2 (f ( 24 ( r2 - 1) ^ 2 + 48?/ (r 2 - 1) n + 48tr 2 



+7r 2 + 24n /2 (r 2 - l) - 9) y 4 + 2fi L t (12 (3r 2 - l) 



n 2 
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+24n' (3r 2 - l) n + 2Atr 2 + 13r 2 + 12n' 2 (3r 2 - l) - 3) y 2 + Q 2 L (3 - 5r 2 ) ) 

+ , f 4 (t 2 (-6 (3r 2 - 2) n 2 + 12nV 2 n + 7r 2 

-6n' 2 (3r 2 - 2) - 9) y 4 - 2VL L t (6 (3r 2 - 2) n' 2 + (l8n 2 - 7) r 2 + 9) y 2 

(6 (5r 2 - 2) n 2 - 12n'r 2 n + 7r 2 - 6n' 2 (3r 2 - 2) - 9) ) ) ) , (C.4) 

(n(t-tt L )+n'(tt L + t))T 

Vt(n L + t) 

(n + n')(2n 2 -An'n + 2n' 2 -At-l)r 2(ra + n')-v/*T 
+ 2~Vt V^i + t 2 ' ' 

(n + n')(n(t - ft L ) + n'{Q L + t)) ( 2n 2 - An'n + 2n' 2 -At- 4f - 1 J r 2 

2t(fi L + t) 

(n + n') 2 (-2n 2 + An'n - 2n' 2 + At + 1) r 2 



- ( (-6 (3r 2 - 2) n 2 - 12n'r 2 n - 5r 2 



i2t(n L + ty 

+6n' 2 (5r 2 - 2) + 3) Q 2 L + 2t (6 (5r 2 - 2) n' 2 - (l8n 2 + 5) r 2 + 3) fi L 
+t 2 (6 (5r 2 - 2) n 2 + 12nV 2 n - 5r 2 + 6n' 2 (5r 2 - 2) + 3) ) 

+-|- ( - 12(n + n') 2 (n 2 - 2n'n + n' 2 - 2t) r 2 

+3 ((5r 2 - 2) n 2 + 2n' (5r 2 - 2) n + At (r 2 - l) + n' 2 (5r 2 - 2)) 
+2 (-7 (3r 2 - 1) n 4 + 4n' (3r 2 - l) n 3 + 6 ((3r 2 - l) n' 2 
+2t (r 2 - 1)) n 2 + 4n' ((3r 2 - l) n' 2 + 6t (r 2 - l)) n 

+Q(n + n') 2 (n 2 - 2n'n + n' 2 - 2t) 2 r 2 + 12t 2 r 2 + 12n' 2 t (r 2 - l) 
-7n' 4 (3r 2 -l))) 

+ ^ f (96t 2 r 2 + 24 (2n 2 + An'n + 2n' 2 + l) t (r 2 - l) 

A8t(u 2 L + t*f 2 \ K V ) \ ) 

+ y/ul + t* (5r 2 -3))cl 

+t 2 (96tV + 24 (2n 2 + An'n + 2n' 2 + l) t (2r 2 - l) 

+ sJuj 2 L + t 2 ((96n 2 + 192n n + 96n /2 + 53) r 2 - 3)^ ^, (C.6) 
2(n + n')-\/*T {uj 2 l + tV) y 2 (n + n') (2n 2 - An'n + 2n' 2 - At - 1) r 



V^i+^K-tV(y 2 -2)) 2Vt 
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r (nty 2 + n'ty 2 - nVt L + n'VL L ) 

Vt(t y 2 + Q L ) ' 1 j 

(i) (n + n') 2 {-2n 2 + An'n - 2n' 2 + At + 1) r 2 (w£ + i 2 y 2 ) y 2 



(e<5)TM 



^K-*V(» 2 -2)) 



n' 2 



+2n' (5r 2 - 2) n + At (r 2 - l) + n' 2 (5r 2 - 2)) 
+2 (-7 (3r 2 - 1) n 4 + 4n' (3r 2 - l) n 3 + 6 ((3r 2 - l) 
+2t (r 2 - 1)) n 2 + An' ((3r 2 - l) n' 2 + Qt (r 2 - l)) n 

+6(n + n') 2 (n 2 - 2n'n + n' 2 - 2t) 2 r 2 + 12t 2 r 2 
+12n /2 t (r 2 -l)-7n /4 (3r 2 - l)) ) 

J- 2 ( 2 (* 2 (- 6 ( 3r2 - 2) n 2 + 12n'r 2 n + 7r 2 



2At (ty 2 + Q L 
_ 6n '2 ^ 3r 2 _ 2 ) _ g) y 4 _ 2fiit ( 6 ^ 3r 2 _ 2 ^ n « 

+ (I8n 2 - 7) r 2 + 9) y 2 + Q 2 L (6 (5r 2 - 2) n 2 - I2n'r 2 n + 7r 2 

-6n' 2 (3r 2 - 2) - 9) 



+ 2t {ty 2 + Q L ) ( n + ^ + n ' ty2 ~ nVtL + 
2n 2 -4nW2n' 2 -4t + , ^M^ffl 1 



+ T^T— -77 ( (("96tV 



v^IT^K-tV(y 2 -2)) 

i 

A8t K + t 2 Y /z (J 2 - t 2 y 2 (y 2 -2)) 
-72 (2n 2 + An'n + 2n' 2 + l) ty 2 + 5^u 2 +t 2 ^j r 

+24 (2n 2 + An'n + 2n' 2 + l) ty 2 - 3^u 2 +t 2 ^J 

+t 2 ( (24 (2n 2 + An'n + 2n' 2 + l) t (y 4 - 9y 2 - 4) y 2 
+9Qt 2 (y 4 - 3y 2 - l) y 2 

+ ^u 2 +t 2 (2 (48n 2 + 96n'n + A8n' 2 + 19) y 4 + 20y 2 + 5)^ r 

+3 (y^ 2 + t 2 (2y 4 - Ay 2 - l) 

-8 (2n 2 + An'n + 2n' 2 + l) ty 2 (y 4 - 3y 2 - l))) u 4 L 
+t 4 y 2 (r 2 (24 (2n 2 + Ann + 2n' 2 + l) t (y 4 - Ay 2 - 12) y 2 

+96t 2 (y 4 -y 2 -3)y 2 + ^ uj 2 l + t 2 (5y 6 + 4 (48n 2 + 96n'n 
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+48n' 2 + 19) y 4 + 10y 2 + 20)) 
-3 (8 (2n 2 + An'n + 2n' 2 + l) t (y 4 - y 2 - 3) y 2 

+ ^2 + t 2 (y 6 -4y 4 + 2y 2 + 4))^l 

+t 6 y 4 (r 2 (48 (2n 2 + 4n'n + 2n /2 + l) t (y 2 - 4) y 2 

+96t 2 (y 2 -2)y 2 + ^uj 2 +t 2 ((96n 2 + 192n'n + 96n' 2 + 53) y 4 
-20y 2 + 20)) - 3 (y 2 - 2) (8 (2n 2 + An'n + 2n' 2 + l) ty 2 

+ y/ul + t* (y 2 -2)^Y (C.8) 



Appendix D 

In this appendix we collect the formulas for the functions f\, introduced in Eq. (!86|) 
describing the contributions beyond PFA. For the considered models these functions 
read 

dt t 3 ' 2 gM™(QL,t) , 

dtt 3 ' 2 I dy gW™(nL,t,y) , 
Jo 

dt t 3/2 g {e5) ™(nL,t) , 

dtt 3 ' 2 I dy g^™(nL,t,y) , 



f[ 5S)TB (QL) 


A80V2 36 
vr9/2 7 


POO 

Jo 


fi 55) ™(ttL) 


480^ 


1 


vr 9 / 2 1- 

36 


40 

37? 


f^(QL,u p L) 


480^ 36 
vr 9 / 2 7 


/■oo 

Jo 


f^™(nL,cu p L) 


480^ 


1 


vr 9 / 2 1- 

36 


40 

3tt 2 



where the functions g for the (<W)-model are given by 
(I6t 4 + 32tt L t 3 + 32t 3 + l6tt 2 L t 2 + 32tt L t 2 + 16t 2 ) Zi_f 



e - 2t n 2 L 



e- 2t Q 2 



+ (-40t 3 - 80Q L t 2 - 24t 2 - 40fi|i - 40Q L t) Li_i f — + ^ j 
+ (32t 4 + 64Q L t 3 + 16t 3 + 32fi 2 t 2 + 16Q L t 2 + lit 2 + 30Q L t + 15fi|) Lii 
+ (-8t 3 - 16Q L t 2 + 12t 2 - 8ft£t + 16Q L t) Lii I 1 

+ (-3Q 2 - 6^ L - 3t 2 ) Lis (^f^J > (D- 1 ) 
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gVQ™(nL,t,y) = 




-8t ( t (-3y 2 + t(y 2 + 3)-l)y 4 

+n L {-hy 2 + 2t (y 2 + 3) + l)y 2 + Q 2 L (y 2 + 3)) Li_i 



(ty 2 + V L ) 2 



+ (-32*V + 16tV - 23t 2 y 6 + 32t 4 y 4 - 64fi L *V - 16t 3 y 4 + l6tt L t 2 y A - A9t 2 y A 
-54tt L ty 4 + 64tt L t 3 y 2 - 27tt 2 L y 2 - 32Vt\t 2 y 2 - l6Q L t 2 y 2 - 90Q L ty 2 



The corresponding function for the (e5)-model have a similar structure but they are 
too lengthy as to be displayed here. 

References 

[1] E. M. Lifshitz. The theory of molecular attractive forces between solids. Soviet 
Physics JETP-USSR, 2(l):73-83, 1956. 

[2] B. Derjaguin. Kolloid Z., 69:155, 1934. 

[3] M. Bordag. The Casimir effect for a sphere and a cylinder in front of plane and 
corrections to the proximity force theorem. Phys. Rev., D73: 125018, 2006. 

[4] H. Gies and K. Klingmueller. 

[5] Aurel Bulgac, Piotr Magierski, and Andreas Wirzba. Scalar Casimir effect 
between Dirichlet spheres or a plate and a sphere. Phys. Rev., D73:025007, 



[6] T. Emig, R. L. Jaffe, M. Kardar, and A. Scardicchio. Casimir interaction 
between a plate and a cylinder. Physical Review Letters, 96(8):080403, 2006. 
cond-mat/0601055. 

[7] M. Bordag. The Casimir effect for thin plasma sheets and the role of the surface 
plasmons. J. Phys. A: Math. Gen., 39:6173-6185, 2006. 

[8] Gabriel Barton and Claudia Eberlein. Plasma spectroscopy proposed for Cqo 
and C 70 . J. Chem. Phys., 95:1512-1517, 1991. 




(D.2) 



2006. 



33 



[9] G. Barton. Casimir effects for a flat plasma sheet: I. Energies. J. Phys., 
A38(13):2997-3019, 2005. 

[10] M Bordag, B. Geyer, G. L. Klimchitskaya, and V. M. Mostepanenko. Lifshitz- 
type formulas for graphene and single-wall carbon nanotubes: van der Waals 
and Casimir interactions. Physical Review B, to appear, 2006. 

[11] M. Bordag, D. Robaschik, and E. Wieczorek. Quantum Field Theoretic Treat- 
ment of the Casimir Effect. Ann. Phys., 165:192, 1985. 

[12] H. Li and M. Kardar. Fluctuation-induced forces between manifolds immersed 
in correlated fluids. Physical Review A, 46(10):6490-6500, 1992. 

[13] T. Emig and R. Buscher. Towards a theory of molecular forces between de- 
formed media. Nuclear Physics B, 696(3) :468-491, 2004. 

[14] O. Kenneth and I. Klich. Opposites attract: A theorem about the Casimir 
force. Physical Review Letters, 97(16), 2006. 

[15] M. Bordag, D. Hennig, and D. Robaschik. Vacuum energy in quantum field 
theory with external potentials concentrated on planes. J. Phys. A, A25:4483, 
1992. 



34 



